NPS ARCHIVE 
1968 

MCMILLAN, J. 






TEMPERATURE FLUCTUATIONS 
NEAR THE AIR-SEA INTERFACE 

by 



John Gabriel McMillan 




library 

MVAL POSTGRADUATE SCHOOfc 
viOlfTEREY. CALIF. 93940 




# 



I 




J 












f 



■s 



I 








UNITED STATES 

NAVAL POSTGRADUATE SCHOOL 




THESIS 

TEMPERATURE FLUCTUATIONS 
NEAR THE AIR-SEA INTERFACE 

by 

John Gabriel McMillan 
June 1968 




TR Citatfoii Database 



Page l-of2. 



Add to Shopping Cart | Save for Later | 
Technical Reports Collection 

Citation Format: Full Citation (IF) 



Accession Number : 

AD0855533 
Citation Status: 

Active 

Citation Classification: 

Unclassified 

SBI Site Holding Symbol: 

AWS/TECH 
Field(s) & Group(s): 

080300 - PHYSICAL AND DYNAMIC OCEANOGRAPHY 
140200 - TEST FACILITIES, EQUIPMENT AND METHODS 
Corporate Author: 

NAVAL POSTGRADUATE SCHOOL MONTEREY CA 

Unclassifled Title: 

Temperature Fluctuations Near the Air-Sea Interface. 

Title Classification: 

Unclassifled 
Descriptive Note: 

Master's thesis, 

Personal Author(s): 

McMillan, John Gabriel 
Report Date: 

01 Jun 1968 
Media Count: 

53 Page(s) 

Cost: 

$7.00 

Report Classification: 

Unclassified 

Descriptors: 

(*SEA WATER, (^METEOROLOGICAL INSTRUMENTS, , SURFACE TEMPERATURE), 
DESIGN), THERMISTORS, CALIBRATION, INTERFACES, WIND, VELOCITY, POWER 
SPECTRA, PERIODIC VARIATIONS, THESES 

Identifiers: 

AIR WATER INTERACTIONS, *AIR SEA INTERFACES. 

Abstract: 

The design, construction and calibration of suitable instrumentation for measuring temperature 
fluctuations within a few centimeters of the sea surface, surface wave variations, and wind 
velocity are discussed. Small bead-in-glass thermistors mounted at 2.5 and 5 centimeters above 
the water surface were used to measure temperature variations under varying wind conditions and 
spectral densities computed on a digital computer. Spectra, co-spectra and coherence squares from 
0 to 5 Hz are presented. (Author) 

Abstract Classification: 

Unclassifled 

Distribution Limitation(s): 



https://drols.dtic.mil/search97cgi/docview.dll?Key=AD0855533&Query=%28%28mcmillan... 4/9/2003 



TR Citation Database 



Page 2-of 2 • 



01 - APPROVED FOR PUBLIC RELEASE 
Source Code: 

251450 

Document Location: 

DTIC 

Change Authority: 

ST-A PER USNPS LTR, 1 MAR 72 



https;//drols.dtic.mil/search97cgi/docview.dll?Key=AD0855533&Query=%28yo28mcmillan... 4/9/2003 



NAV^POSTGRADUATE schooU 
MONTEREY. CALIF. 93940 



TEMPERATURE FLUCTUATIONS 



NEAR THE AIR- SEA INTERFACE 



by 



John Gabriel McMillan 
Lieutenant Commander, Ufiited States Navy 
B, S., Naval Academy, 1958 



Submitted in partial 
requirements for 



fulfillment of 
the degree of 



the 



MASTER OF SCIENCE IN OCEANOGRAPHY 
from the 

NAVAL POSTGRADUATE SCHOOL 



June 1968 



ABSTRACT 



The design, construction and calibration of suitable 
instrumentation for measuring temperature fluctuations 
within a few centim.eters of the sea surface, surface wave 
variations, and wind velocity are discussed. Small bead- 
in-glass thermistors miounted at 2.5 and 5 centimeters 
above the water surface were used to measure temperature 
variations under varying wind conditions and spectra] 
densities comiputed on a digital computer. Spectra, co- 
spectra and coherence squares from 0 to 5 Hz are presented. 
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Introduction 



The energy exchanges between the atrr.o sphere and the 
oceans are the driving force for a therrrodynaric and rrech- 
anical systeir which enconpasses the fluid portions of the 
earth. The atmosphere derives energy from the sea in the 
form of sensible and latent heat and in turn, provides 
mechanical energy for oceanic circulations. Sm.all-scaie 
processes, such as convection and local wind fields, are 
also direct results of these energy exchanges, and the 
ability to accurately measure energy fluxes at any level 
over the sea would permit one to predict the growth and 
decay of such processes. 

There are generally four forms of air-sea energy 
exchange that one may consider: 

1. Electrom.agnetic radiation 

2. Convection and conduction of sensible heat 

3. Convection and diffusion of latent heat in the 
for, of water vapor 

4. Vertical exchange of mechanical energy. 

In investigating the flux of any of these forms of energy, 
one n.ust measure directly certain quantities. For heat 
flux, the appropriate equations are; =^Cpfw and E^ = 

/O Lw, where E is the sensible heat flux. En is the latent 
heat flux,^and^^^ are the densities of air and water re- 
spectively, Cp is the specific heat of air at constant 
pressure, L is the latent heat of evaporation (or conden- 



sation), T is the tempt ratlin and w is the vertical cot: - 
ponent of the wind velocity. Then, to measure the heat 



fluxes, one roust be able to Treasure tenperature , vapor 
pressure, and wind velocity at any desired level. Since 
waves on the water surface interact with the wind field 
(which carries the scalar temperature field), the sur= 
face wave field itay also have an effect on heat flux. 

As a preliminary step in the determination of heat 
flux, methods of m-easuring variations of temperature, 
wind velocity, and water surface displacement on the 
smallest possible time and space scales m.ust be devised 
and tested. This study was made to investigate some of 
the possible approaches to these m.easuren;ents close to 
the air-sea interface and had as its objectives: 

1. The design of a suitable sensor system for 
measuring the temperature fluctuations within 
a few centimeters of the sea surface, the wave 
field at the same point, and the downstreamj 
wind velocity at a suitable level above the 
surface . 

2. The continuous recording of these three meas- 
urer, ents in a form suitable for digital anal- 
ysis, 

3. The analysis of the records to determine the 
appropriate spectra and the relation, if any, 
between the air- temperature variations and the 
water- surf ace height under varying wind con- 
ditions . 
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Design of the Measu re -rent. S yste n. 

A systeij to n'easure srrail temperature variations 
requires a temperature sensor that is durable, inexpen- 
sive, and that has a sufficiently small response time to 
record fluctuations at frequencies of about 10 Hz. A 
small bead thermistor was used and although thermistor 
time constants are generally large compared to, say, a 
platinumj resistance thermometer, it appeared that they 
were sufficiently small for this investigation. A Vic- 
tory Engineering Corporation bead-in-glass thermistor, 
model 32A49 , was used. The manufacturer's specifica- 
tions for this unit are. a nominal diameter of 0.013", 
a time constant ('T') of 1.0 seconds in still air, a pow- 
er dissipation constant of 0.10 mw/'^C, and a continuous 
operational range of -50°C to 300'^C. The zero-power 
resistance at 25°C is 2500XL ^ 25%, and the resistance- 
tern, perature ratios are 5.6 for at O'C/R^ at 50°C, and 
14.6 for R^ at 25°C/R^ at 125 °C. However, over the range 
of temperatures encountered in this study f about 5°C), the 
resistance-temperature curve is essentially linear (Figure 

7). 

Two of these therrristors were mounted on a guided 
styrofoam float, at about 2.5 and 5 cm above the water 
surface, as shown in Figure 2. A variety of float designs 
were tested for response to surface waves. The one used 
in the field was a hemisphere-" of about 2.5 err diameter. 

In wave tank tests, movies showed the float to follow the 
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waves quite well, although, of course, it tended to act 
as a filter for most waves smaller than the diameter of 
the float. 

To measure the water-surface fluctuations, a resist- 
ance wave gauge was mounted as close as possible to the 
thermistors. This gauge consists of a pair of 0,01” di- 
ameter stainless-steel wires, about 25 cm in length and 
positioned some 3 cm apart (Figure 1). A 5 volt, 20 Hz 
power source was used to avoid electrolysis, and the sig- 
nal output from the bridge was fed through a full-wave 
rectifier . 

Wind velocity was measured with a Caseila sensi- 
tive cup anem.ometer, mounted approximately one-half meter 
above the water. The anemometer has a digital counter to 
record wind velocity and was used to determine mean wind 
speeds over time intervals of 10 to 30 seconds. 

The wave gauge and thermistors were installed as 
one leg of separate Wheatstone bridges and the outputs 
recorded on Mosley 7100B strip recorders. This recorder 
has individual DC amplifiers for each channel with maxi- 
mumj amplification providing a 5-mv per 10-inch scale. 

The recorders were operated at 0.5 inches-per- second chart 
speed for all measurements. The bridge and rectifier cir- 
cuitry for each instrument is shown in Figures 3 and 4, 
and a block diagram of the basic data-collection system 
is given in Figure 5. 
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The platfori:, for fit id ounLiog Lht sensors is 
shown installed for neasurer ents in Figure 6, The trans 
mission lines for signal output are shielded, two-core 
telemetry cables with low-loss characteristics. 

^ * Instrument Calibration a nd Observation Tec hniques 
The therm, istors were calibrated by varying the tern 
perature of a water bath and plotting the resulting re- 
sistance change for each temperature increment. The re- 
sults a.re presented in Figure 8. Since this curve was 
essentially the sane for t ach indi\ idual thermistor, it 
was assumed that minor differences in slope could be ig- 
nored. Before field measurements, the bridges were bal- 
anced in the laboratory at a temperature close to that 
expected in the field. Each recorder channel was cali- 
brated to a scale of 1°C per inch and for equal sensi- 
tivity in each channel. The entire system, was then test 
ed with a known temperature change with the telemetry 
cables installed. The response time of the system, was 
estimated by running the recorder at its maximum, chart 
speed (2 inches per second) and recording a known temper 
ature change. The time for 63% of the change (one time 
constant) was then measured directly from the record. 
Since the temperature changes were induced by adding 
water of a different temperature to the bath, the actual 
temperature change was not instantaneous, and the esti- 
mate of 80 msec is perhaps too large. It is, however, a 
good indication that a flat response to 10 H? is a con- 
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servative estimate. This response will be somewhat less 
in mioving air, and a response of 5 Hz was considered to 
be the upper limiit for a flat response for this study. 

The wave gauge was calibrated in a wave tank to find 
the required immersion for maximum signal amiplitude. The 
smaller the immersion, the greater the signal amplitude. 
Above a wave amplitude of about 5 cm, the output became 
slightly distorted due to the non-linearity of the gauge 
system. Again, the range of measurements in this study 
was sufficiently sm.all that the response could be con- 
sidered linear. The output of the full-wave rectifier 
was compared with a comm, on generated signal of 0 to 10 Hz 
and showed an undistorted output over this range. 

The Case 11a anemor. tier vvas not calibrated. The manu- 
facturer's conversion curves for counter readings to veloc- 
ities were assumed to be accurate. 

The field observations were made during the months 
of April and May, 1968 at Roberts Lake, Seaside, Califor- 
nia. This small lake is approximately one-half mile in- 
land from Monterey Bay and is about 450 yards long and 
200 yards wide. Bctvceen the bay and the lake are a series 
of sand dunes about 25 to 50 feet high which end about 150 
yards from the lake. There is no other obstruction to air 
flow in this area to seaward and the winds are essentially 
those present over Monterey Bay. The available unrestrict- 
ed fetch at the measurement site is between 75 and 400 
yards, depending on the wind direction. Most of the data 
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was taken with a fetch of frora 100 to 300 yards. The 
platform was placed in about three feet of water, 25 to 
30 feet from shore, and the signals were trans'^^itted to 
the recording equipment with the low- loss cables. 

Measurements were made on six occasions, the first 
three with only a single thermistor mounted at 2,5 cm. 

The lengths of the recordings were between 3 and 5 min= 
utes, with mean measured winds from 2.4 to 8 meters per 
second. The significant wave heights were between 2 
and 4 cm, with periods of from 1.2 to J, .8 seconds and 
lengths of fromi 10 to 15 cm. One representative record 
of the single-probe measurements and one of the dual- 
probe measurements were analyzed. A portion of a compu- 
ter di'awn record of dual -probe data is shown in Figure 8. 

^ • Ana lysis 

In previous air- temperature measurements (Pond et 
al , 1966), a substantial portion of the spectral "energy” 
was concentrated in the region below 10 Hz, Since 5 Hz 
was considered to be near the upper limit of the flat 
response of the thermistor, the analysis of the data was 
restricted to the range of frequencies from 0 to 5 Hz, 

As the data was recorded in a graphic form, the conver- 
sion to digital form required a manual interpretation at 
intervals of 0,05" and 0,025'^C, This process certainly 
introduced some noise into the record (,in addition, to 
that generated by the electronics in the system), but this 
was unavoidable. The ensuing analysis showed this noise 
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to be negligible. 

Since the tertperature fluctuations at the air-sea 
interface are felt to be essentially randon and station- 
ary, it is appropriate to analyze the data by computing 
power spectra. This irethod of analysis is the basis of 
an existing FORTRAN IV prograrrj in the library of the 
N.PoGoS. computer library, F^ROGRAM BLACKY (Table I) 

This program corr-putes smoothed power spectral estiirates 
frorr the Fourier transforrr of the covariance for a given 
frequency bandwidth. 

In analyzing any continuous process using discrete 
points, the problerr: of aliasing is encountered. For this 
study, it was assurr.ed that little energy existed above 5 
Hz. PROGRAM BLACKY also provided a filter by running a 
weighted average for thre^e consecutive points of the 
record. This would reduce the possibility of folding and 
also reirove most of the high frequency noise.- 

In order to confirm the r« suits of PROGRAM BLACKY, 
and to obtain a graphic display of the spectra, the data 
was also run in PROGRAM BIMED 02T, an autocovariance and 
power spectrum program of the Health Science Facility, 
UCLA The records were also analyzed for means, density 
and cum^ulative distributions. 

5 . Results and Conclusions 

The data was processed for smoothed estimates of 
spectral density from. 0 to 5 Hz at an interval of 0,025 
cycles. The record was filtered by means of a Hanning 
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lag window of the forrr, : 



D =0,5x1+ cos 
r 



where is the weighting function, id the nuTi'ber of lags 
for the record, and r is an increiDental index of m . The 
weights for this analysis were 0.25, 0,50 for the three 
consecutive points. The first and last points were com- 
puted with weights of 0.5 + 0.5 x^^ and 0.5 + 0,5 

x^. This low-pass filter effectively removed the high 
frequency noise of the record. Since little energy was 
felt to exist at high frequencies, the loss of this por- 
tion of the signal was considered to be insignificant. 

The smoothed spectral estimates for the temperature 
variations are presented in Figures 9 and 10, and the 
corresponding wave amplitude spectra in Figure 11. A co- 
spectrum for the dual-probe temperatures is presented in 
Figure 12. The values from which these plots were taken 
are given in Tables II and III. The spectra are presented 
as log^Q^^xf) versus logj^Q(f), with *+^(1) in units of 
(°C) sec. The spectrum is defined such that 



where 9 is the temiperature variation iBatchelor, 1959). 
The wave amjplitude spectrum, 0(f), is defined similarly. 



zen turbulence" vTaylor, 1938), is not v^alid for a region 
of shear flow, the relation 




2 

with units of cm sec. Since Taylor's hypothesis of "fro- 



k 
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jL / « 



2JLJ 

u 



where k is the wave nurrbcr, i the frequency, and U the 
mean wind velocity, was not used and the spectra ar^ pre- 
sented in terms of frequency rather than wave number,. 

The total "energy” was estimated in two ways in PRO- 
GRAM BIMED 02T; (1) by computing the variance, and (2) by 
integrating the power speccrum with the trapezoidal rule. 
The ratio of (l)/(2; giv'cs the traction of the total ener- 
gy in the computed spectra. The following table gives the 
individual values and the resulting ratios. 



t(f)/0(£) 


Check sum 
of PSD 


Energy by 
Trap, rule 


Dif f . 


Ratio 

(i)/(2)_ 


Temp, 2.5 cm 
(Single-probe 


1.1057844 


1.1058750 


-9.06x10^^ 


] .00 


Temp. 2,5 cm 
(Dual-probe) 


0.23003.38 


0.2300585 


2.47x10*^ 


0.99 


Temp. 5 cm 
(Dual-probe ) 


0.3806309 


0.3806635 


3.27x10'^ 


0.99 


Wave Amplitude 
(Single probe) 


2.5749149 


2,5751486 


2.337x10“^ 


0.99 


Wave Amplitudf 
(Dual-probe ) 


0.2900739 


2.2900738 


1,00x10“^ 


1 ,00 



The high percentage of energy in the computed spectra indi- 
cates that there is little f.noigy above 5 Hz and there was 
negligible energy folding from higher frequencies. 

Confidence limits for the spectral estimates were 
calculated for an 80% interval , For the analyzed record 
length of 120 seconds, a sampling interval of 0,1 seconds, 
and a correlation lag value of 200, the 80% limits are; 



0,65 ^(f)^ TCfJ ^ The coeputaLion of these 

limits is presented in Appendix A, 

Coupling, between the surface waves and the tern'" 
perature fluctuations at fixed heights above the undu- 
lating surface, does not appear to exist, even at the 
higher wind speeds. The coherence square function is 



(Bendat and Piersol, 1966). The coherence square for the 
temperature and wave spectra (Table IT) shows virtually no 
correlation between these qiaantlties. Since the tempera- 
ture measurements were taken at a height that remained 
constant with respect to the water surface, this lack of 
correlation suggests that the temperature structure at 
this level is not affected by the passing wave surface. 

In contrast, spectral densities for temperature measure- 
ments taken at a fixed height above mean sea level have 
shown periodic maximja at the same frequencies as wave 
amplitude spectral maxima (Makova, 1963). A comparison 
of m.easure::ents taken at a fixed height above the surface 
and simiultaneously at a fixed height above mean sea level 
would be useful in further studies of this problem. 

The one -dimensional spectra of temperature fluctua- 
tions plotted on a log-log scale (Figures 9 and 10) show 
a / 2 slope for frequencies above 0.2 Hz, This slope has 



f (f) = 



2 




xy 



G if) G vf) 
X y ^ 
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also been found for teir>perature ~casureT:ent s taken with a 
platinum resistance thermometer at heights of one meter 
and greater (Pond et al, 1966). Since the Inequality, 
kz»4.5 does not hold near 0.2 Hz, this result indicates 
that a relation of the form 

y ' f ) - CjAe^e ^ 3 k ^ ^3 

(Batchelor, 1959; is valid t^r ancvalously large scales, 
(Here, is a conversion constant to relate the spretrur 
to frequency rather than wave number, Cq is the scalar 
dissipation, e is the energy dissipation, and A is a uni- 
versal constant.) This result was also found by Pond et 
al (1966). The high values of the water-wave spectra be- 
low 0.05 Hz are not readily explainable, but may be re- 
lated to a periodic pile-up of water at the leeward end 
of the fetch, where the measurements were taken. In a 
small lake, such as was used in this study, periodic 
gusts can easily set up a seiche. 

An exa"'ination cf the coherence square for the dual- 
probe temperatures shows a high degree of correlation over 
the entire spectrum (Table III;, This, together with the 
simularity in the amplitudes of the fluctuations as seen 
in Figure 8, is som.ewhat surprising in view of the 5-cm! 
probe being twice as far from the water as the lower 
probe. The mean temperatures tor thf dual probe record 
are; 14,26 C at 2.5 cm and lu,C4'C at 5 cm. The mean 
air tem:perature as recorded at approximately 5 cm above 
the water with a m.ercury thermometer was 14.0°C and the 
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water terrperaturt at the surface u-as 2Q.r C. Clearly, 
the air near the water was highly unstable, and it ap- 
pears that the greatest air temperature gradient occurs 
be low 2.5 cm . 

The results of this study indicate that the data 
collection system used for obtaining temperatures and 
wave records was reasonably accurate and chat additional 
work should be done to rrduce the lev^el at which the tem- 
peratures can be taken. 

^ • Re co mrrendations for Futu re. S tudi es 

Further work in this field should include a con- 
tinuous record of the wind field, additional thermistor 
probes at lower levels (and at least one below the water 
surface), a simultaneous record of temperature at a fixed 
height abov'e mean water level, and the development of a 
sensor mounting that will not interfert with the wave 
field, such as a servo - fol lower syster;-. 
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APPENDIX A 



Confidence I nterval Calculations 
Record length = = 120 seconds 

Sampling interval = f = 0,1 seconds 
No, of lags per record length = rr = 200 
Degrees of freedom = d 
No. of data points = N 



120 
^ (?:i 



= 1200 



d 




2A00 



= 12 



(Blackman and Tukey, 

1958) 



The confidence interval is defined 



^ /-2 ^ d G~^ 

(Bendat and Piersol, 
Ad,«^/2 Ad,l-®‘/2 1966) 



Assuming that the Fourier transform bears the same reia- 

2 

tion, o< = 0.2 and d = 12, gi'*cs values for Xd*«^/2 

* Xd ; l -«</2 

7d;<=</2 = 26,22 
X^,l.C/2 = 6.30 



^2 
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